Khan and Akbar SpringerPlus 2014, 3:324 
http://www.springerplus.eom/content/3/1/324 



0 SpringerPlus 

a SpringerOpen Journal 



RESEARCH 



Open Access 



Study of analytical method to seek for exact 
solutions of variant Boussinesq equations 



Kamruzzaman Khan and M Ali Akbar 



* Correspondence: 

k.khanru@gmail.com 

'Department of Mathematics, Pabna 

University of Science and 

Technology, Pabna 6600, 

Bangladesh 

Full list of author information is 
available at the end of the article 



Abstract 

In this paper, we have been acquired the soliton solutions of the Variant Boussinesq 
equations. Primarily, we have used the enhanced (G'/G)-expansion method to 
find exact solutions of Variant Boussinesq equations. Then, we attain some exact 
solutions including soliton solutions, hyperbolic and trigonometric function 
solutions of this equation. 
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Background 

The theory of nonlinear evolution equations (NLEEs) has become one of the most 
important fields of study in mathematical physics and engineering. This is essentially 
due to the frequent occurrence of NLEEs in many branches of physics, mathematics, 
and other branches of sciences. With much interest and great demand for applications 
of NLEEs in various areas of sciences, several analytical methods to find exact solutions 
of NLEEs have been developed by diverse group of mathematician and physicist. It is 
significant that many equations of physics, chemistry, and biology contain empirical 
parameters or empirical functions. Exact solutions allow researchers to design and 
run experiments, by creating appropriate natural conditions, to determine these 
parameters or functions. A variety of powerful methods, such as the sine-cosine 
method (Wazwaz 2004a; Bekir 2008), The first integral method (Bekir et al. 2014; 
Jafari et al. 2012), the homotopy perturbation method (Mohyud-Din and Noor 2009; 
Mohyud-Din et al. 2011), the (G VG)-expansion method (Wang et al. 2008; Zayed 
and Gepreel 2009; Guo and Zhou 2010), the Exp-function method (Bekir and Boz 
2008; Akbar and Ali 2011; Naher et al. 2010; Ebadi et al. 2013), the modified simple 
equation method (Jawad et al. 2010; Zayed 2011; Khan and Akbar 2013a, b), the exp 
(-<2>(£))-expansion method (Khan and Akbar 2013c), the Enhanced (GYG)-expansion 
Method (Khan and Akbar 2013d, e; Islam et al. 2014), the tanh-function method 
(Wazwaz 2004b, 2005, 2007), and the modified tanh-coth function method (Lee and 
Sakthivel 2011) were used to find new exact traveling wave solutions of nonlinear 
evolution equations in mathematical physics. 

Various ansatz have been proposed for seeking traveling wave solutions of nonlinear 
differential equations. Recently, Wang et al. (2008) have introduced a simple method 
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which is called the (G7G) -expansion method to look for traveling wave solutions of 
nonlinear evolution equations, where G = G(£) satisfies the second order linear ordin- 
ary differential equation G" +AG' +pG = 0, where A and [i are arbitrary constants and 



= a m \^rj + be the traveling wave solution of NLEEs. By means of this 

method they have solved the KdV equation, the mKdV equation, the variant Boussinesq 
equations and the Hirota-Satsuma equations. Guo and Zhou (2010) have intro- 
duced an another method so called extended (GVG)-expansion method where 
G = G(£) satisfies the second order linear ordinary differential equation G" + ^G = 0, 

where G = dG ^p , G = d , £ = x-Vt, V is a constant and u(Q = ao + 
Vj «; (G /G)' + b t (G /G 



i=l 




be the traveling wave solution. 



They proposed extended (G7G) -expansion method to construct travelling wave 
solutions of Whitham-Broer-Kaup-Like equations and coupled Hirota-Satsuma 
KdV equations. 

Among those approaches, an enhanced (G7G) -expansion method is a powerful tool 
to reveal more general solitons and periodic wave solutions of NLEEs in mathematical 
physics and engineering. The focal ideas of the enhanced (G /G) -expansion method 
are that the traveling wave solutions of NLEEs can be expressed as rational functions of 
(G /G), where G = G(£) satisfies the second order linear ordinary differential equation 
G" + \a G = 0. 

The objective of this article is to present an enhanced (G /G) -expansion method to 
construct the exact solitary wave solutions for NLEEs in mathematical physics via the 
Variant Boussinesq equations. 

The article is arranged as follows: Methodology, Application of the enhanced (G7G)- 
expansion method, Graphical representation, Comparisons and conclusions. 

Methodology 

In this section, we discuss an analytical method, so called enhanced (G'/G) -expansion 
method, for deriving traveling wave solutions to PDEs. We will first discuss the method 
applied to a problem defined in terms of a nonlinear partial differential equation having 
two independent variables, one space dimension x, and another the time dimension 
t. Subsequently, it will be shown that the arguments extend naturally to coupled 
equations and also to problems defined in terms of two or more spatial dimensions, 
plus time. 

Suppose the evolutionary equation, say in two independent variables X and t, for 
which we wish to find traveling wave solutions is given by 

u t , u x , u t t, u xx , u xt , ) = 0, (1) 



where w(£) = u(x, t) is an unknown function, 91 is a polynomial of u(x, t) and its partial 
derivatives in which the highest order derivatives and nonlinear terms are involved. In 
the following, we give the main steps of this method (Khan and Akbar 2013a, 2013b; 
Islam et al. 2014): 
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Step 1. Combining the independent variables X and t into one variable £, we suppose 
that 

u(Q = u(x,t),Z = x±Vt 1 (2) 

where Ve 91 - {0} is the wave velocity. 

The traveling wave transformation Eq. (2) permits us to reduce Eq. (1) to the follow- 
ing ODE: 



p{u,u ,u , ) =0, (3) 

du ' ' ( c\ d 2 u 

dv u de 



where p is a polynomial in u(£) and its derivatives, while u (£) = % , u (£) = ^3 and 



so on. 

Step 2. We suppose that Eq. (3) has the formal solution 



= - / a+G'/G)' b , G - /G y-i 




(4) 



where G = G(£) satisfy the equation 

G +^G = 0, (5) 

in which a b b t (-« < i < n; n e N) and A are constants to be determined later, and a=± 
\,H * 0. 

Step 3. The positive integer n can be determined by considering the homogeneous 
balance between the highest order derivatives and the nonlinear terms appearing in 
Eq. (1) or Eq. (3). Moreover precisely, we define the degree of w(<f) as D(u(£)) = n which 
gives rise to the degree of other expression as follows: 

D {^)= n+q > D { u i^))= np+s{n+q) (6) 

Therefore we can find the value of n in Eq. (4), using Eq. (6). 

Step 4. We substitute Eq. (4) into Eq. (3) using Eq. (5) and then collect all terms of 



same powers of (G7G/ and (G /G)'J a {\ + i (G /G) 2 ^ together, then set each coeffi- 
cient of them to zero to yield an over-determined system of algebraic equations, solve 
this system for a it b b A and V. 

Step 5. From the general solution of Eq. (5), we get 

When fi < 0, 

— = ^ tanh^ + y/=i$) (7) 

And — = y^cothU + v^f) (8) 
G 

Again, when ft > 0, 

|-=V^tan(A-^) (9) 

And — =y/ficot(A+^) (10) 
G 
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where A is an arbitrary constant. Finally, substituting bi (-« <i<n;n& N), A, V and 
Eqs. (7)-(10) into Eq. (4) we obtain traveling wave solutions of Eq. (1). 

Application 

In this section we will exert the enhanced (GVG)-expansion method to derive the exact 
traveling wave solutions of the Variant Boussinesq equations in the form (Wang et al. 
2008) 

H t +(Hu) x + u xxx = Q, (11) 
u t + H x + uu x = 0, (12) 

where u = u(x, t) represents the velocity and H = H(x, t) represents the total depth of 
water waves. 

The traveling wave transformation equation u = u{x, t), H(x,t), f; = x- Vt converted 
Eq. (11) and Eq. (12) to the following ordinary differential equations for u = u{£) and 
H = H(0. 

-VH +{Hu) +u"' =0 (13) 
-Vu + H + uu = 0 (14) 
Integrating Eqs. (13) and (14) once with respect to <f, we obtain 
R-VH + Hu + u" =0, (15) 

S-Vu + H + ^u 2 = 0, (16) 

where 7? and 5 are constants of integration. 
Now Eq. (16) yields, 

H=Vu-^u 2 -S (17) 

Substituting Eq. (17) into Eq. (15), we get 

u "-(V 2 + S)u + ^Vu 2 -^u 3 + VS + R = Q (18) 

The homogeneous balance between u" and u 3 yields n = 1. 
Hence for n = 1, Eq. (4) reduces to 



g_i(l + X(G'/G)) ai (G/G) J^ + 1 (r■/r^ 2 

+ bo(G' /G) _ y<r(l + X - (G'/G) 2 ) + b^a{\ + i (G /G) 2 ) , 

(19) 

where G = G(0 satisfies Eq. (5). Substituting Eq. (4) along with Eq. (5) into Eq. (3), we 



get a polynomial of (G'/G) 7 and (G /G)' Ja{\ + ^ (G /G) 2 ^ . From these polynomials 



we equate the coefficients of (G'/Gy and (G /G)y ff^l + ^ (G /G) j , and setting 
them to zero, we get an over-determined system that consists of twenty-five algebraic 
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equations. Solving this system for a b bi, A and V, we obtain the following values with 
the aid of symbolic computer software Maple 13. 



Case 1:R = 


0,5 


= 9,pi + \a 0 2 , V 


— dQ,X 


= 0, «_i 


= — 2pi, do = do, d\ = 2, 


= b 0 = 


fcl = 


0, 


Case 2: R = 


0,5 


= 9,pi + \a 2 ,V = 


- do , X 


= 0, a_i 


= 2pi, ao = ao, a i — ~2, b-i 


= b 0 = 


bi = 


0. 


Case 3: R = 


0,5 


= \pi + \a 0 2 , V = 


- CLq.X = 


= 0, «_i = 


- 0, «o = fl 0i fl i = 1, = 


= 0,^i 


= ± V 


/I 

a' 


Case 4: R = 


0,5: 


= \pi + \ao 2 , V = 


Clo,X = 


0,a_i = 


0, #o = ; &i = ~ 1 ) b-i = bo 


= 0,&i 


= ± V 


k 
a- 



Case 5: 7? = 0, S = 2pi 2 X 2 + 2/4 + 2Xpido + |«o 2 , V = 2piX + do, A = X,d-\ = 2ft, 
a 0 = a 0 ,ai = 0,b-i = b 0 = h = 0, 

Case 6:R= -8^ 2 A(l + piX 2 ) , S = -^pi 2 X 2 -^pi + 2Xpido + \do 2 , V = 2piX + d 0 ,X = A, 
d-i = 2/i, da = do, d\ = 2 + 2^A 2 , b-\ = bo = b\ = 0, 

Case 7: R = 0,S = 2pi 2 X 2 + 2pi-2Xpido + \do 2 , V = -2piX + do,X = X, 
d-i = -2pi, do = do, d\ = 0, b-i = bo = b\ = 0, 

Case 8:R = 8pi 2 X[l + pX 2 ),S = -4pi 2 X 2 -4pi-2Xpido + \ao 2 , V = -2piX + d 0 ,X = X, 

d-i = -2pi, do = do, d\ = -2-2piX 2 , b-i = bo = b\ = 0, 
Case 9: 7? = 0, S = -pi + \do 2 , V = do, X = A, «_i = 0, do = do, d\ = 0, b-\ = bo = 0, 
*i = ±2^, 

Case 10: R = 0,S = 2pi 2 X 2 + 2pi + 2Xpido + \do 2 , V = 2piX + do,X = X,d-\ = 0, 

do = do, d\ = 2 + 2pi A 2 , b-i = bo = b\ = 0, 

Case 11: R = 0,5 = 2^ 2 A 2 + 2pi-2Xpido + \do 2 , V = -2piX + do,X = X,d-\ = 0, 

do = do,d\ = -2-2^/A 2 ,^_i = &o = ^i = 0, 

Case 12: R = 0, 5 = + -flo 2 s V = a 0 , A = A, a_i = 0, «o = flo, «i = b-\ = 0, b 0 = ±^ , b\ = 0. 
Case 13: R = 0,5 = -Xpido + \do 2 + \[* 2 X 2 + V = -Xft + do, X = X, a_i = -pi, do = do, 

d X = 0,&_i = Q,bo = ±-h=,h = 0. 

Case 14: R = 0, S = Xpao + \do 2 + \n 2 X 2 + V = Xpi + do, X = X, d-\ = pi, do = do, 
d x = 0,&_i = 0,&o = ±-^=,bi = 0. 

Now putting the values of d it b it X and V into Eq. (19) and Eq. (17), we obtain the 
following exact solutions. 
Hyperbolic Function Solutions for pi < 0: 

Family 1: = fl 0 +2y=^( tanh(^l + + coth(A + y/^pl^)), 

7fi(0 = 2^(tanh 2 (^l + y=££) + coth 2 (A + v^f)-2), 

where £ = x- d 0 t. 
Family 2: w 2 (<f) = d 0 -2^pl( Unh(A + ^plZ) + coth(A + ^pi^)), 
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H 2 tf) = 2^(tanh 2 (,4 + y^<f) + coth 2 (A + y/=j%)-2), 
where f = x - a 0 t. 

Family 3: w 3 (f) = a 0 +^Jl(tanh{A + yf^jiC) + Isech(A + y^f)), 

w 4 (£) = a 0 + v / = ^(coth(A + s/^C) + csc/z(A + y/^f)), 

H 3 (Z) = -fi( sech 2 (A + y/=itf) ±/tanh(A + y/=j%) sech(A + y/=jtf)), 

Httf) = -n{ csch 2 (A + y/=j#) ± coth(A + y/=j^) csch(A + y/=fty), 

where £ = * - a 0 & 
Family 4: «s(£) = fl 0 -y^( tanh(A + <J-Jis\ + /sec/zL4 + y/^ftf;)), 

u 6 {Q = a 0 -y/=fi(coth(A+ y/=fiZ) + csch(A + y^f)) , 

H s (0 = -fi{ sech 2 (A + y/^) + /tanh(A + y^f) sech(A + y/=j%)). 

7f 6 (0 = -p(csch 2 (A + y/^) + Icoth(A + y/=j%) csch(A + y/=f%)), 

where f; = x- a 0 t. 
Family 5: u 7 (Q = a 0 -2 A / = ^( coth(A + y/^/2£) + ^y/-Ji) , 

m 8 (0 = a 0 -2y/^/2( tanh(A + yf^C) + Ay/^) , 

H 7 {£) = 2fi csch 2 (A + y/=j%) , 

H 8 (0 = -2psech 2 (A + ^jiZ)., 

where £ = x - {2fi A + a 0 ) t. 

(2+2^A 2 ) v ^7tanh(A+ v ^f) 2^(l+A v /=^tanh(/t+ ) ) 
Family 6: ) - « 0 + 1+Av ^ tanh ( yl+v ^ f ) + ^tanhl^+^f) ' 

(2 + 2^A 2 )y=^coth(A + yf^jiCj 2^(1 + Ay=^Zcoth(>l + y^O) 
Wio(t) - «o + 1+A y=^ coth ( j4+ yz^c) + y^rath^4Ty=K) 

/ffn = m Y (2 + 2^A 2 ) V ^tanh(A + v ^f) 2^(1 + Ay^tanh(A + y^ff)) 

^ fl o^«o i + A % /r^tanh(A+ y^f) y^tanh(A + 

1 / (2 + 2^A 2 ) % /^Mtanh(A + v / = ^f) 2^(l + Ay^i tanh(A + yf^)) 



2 I 1 + Ay^tanh(A + y^f) y/=^tanh(A + y^f) 

4fi 2 \ 2 + 4^-2A^«o--'5 

, (2 + 2^A 2 ) v ^Scoth(A + y/=tf) 2^(l + A v ^Scoth(A + y^f)) 

"loKJ = 12A^ + «oJ a 0 H : , — = TT^ 1 ^= TT^ 

1 1 + A v / =^coth(-A + y/=flQ y^coth(A + ypjtf,) 



11 (2 + 2^A ) 1 /=^coth(^l+ y^<f) 2^(1 + Ay=^coth(A + yTpCj) 
&o H — - — — — . h 



2 I 1 + Ay^coth(A + y^t/f) v /=^coth(A + y^f) 



2 2 1 ' 

4^ A + 4^-2A^flo _ 2 a o' 



where £ = x - (2A^ + a 0 ) fc 
Family 7: «n(£) = a 0 +2 v /= ^( coth(A + yf-p-C) + Ay 73 ^) 

"12(0 = a 0 +2 v /z ^(tanh( J 4 + y^f) + Ay/^), 
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H 12 (Q = 2f,sech 2 (A+^), 
where £ = x - {do - 1\\a) t. 

_ (2+2 f <A 2 ) v ^tanh(/l+ v ^f) 2p(l+A v ^tanh(A+y=pf)) 

Family 8: u^) - a 0 1+Ay ^ tanh ( A+ ^ f ) y^tanh^+^f) ' 

(2 + 2^A 2 ) 1 /=^coth(A + y^f) 2^(1 + Ay^wcoth(A + v^^)) 
Mi 4 (t) - a 0 i + Ay^coth^ + y=^f) y=^coth(^ + y^^) ' 

(2 + 2^A 2 ) v /=7<tanh(A + v /=^£) 2<w(l + Ay^tanh^A + V^O) 



" " <r 21/ '''l fl ° l+Ay^tanh^ + y^Sf) v C ^ tanh (^ + v^) 

2 + 2/^A 2 )y=^tanh(A + y=^f) 2/*(l + Ay=^tanh(A + V^O) 
2 1 + Ay=^ tanh(A + ^f^E) ^fl tanh(A + y^f) 

4^ 2 A 2 + 4^< + 2A^a 0 - ^«o 2 

(2 + 2^A 2 )y^coth(A + y^f) 2^(l + Ay^coth(A + y^f)) 

"14(C) = («o-2A^) « 0 



' 2 I a ° 1 + A coth (A + y^f) ^ coth(A + y^f) 



1 +Ay=^coth(A + y^f) ^<coth(A + 

(2 + 2^A 2 ) y^coth (A + y^ - ) 2^(1 + A/^coth(A + V^)A 



4^< 2 A 2 + 4^ + 2A^flo _ -«o' 



where £ = x- {a 0 - 1X\a) t. 
Family 9: Wis(<f) = a 0 +2^/fisech(A + y^f), 

w 16 (f) = a 0 ±/2y7<csc/2(y4 + y^f), 

ffis(0 = ^(l-2sec/* 2 (A + y=^f)), 

ifi6(0 = p(l + 2 csc/* 2 (A + v^f) ) , 

where £,=%- a 0 t. 



/ tanh('A+,/=«f) \ 
Family 10, «„(0 = «„+2 v ^(l + ( i^b^fe) J < 



/ sec/2 (A + y^f) 
i 1 + Ay^tanhf/L + y^K), 

csch(A + y/^) X 
l+Ay=^coth(A + v^f) i 



Hi 7 (0 = -2^(1 + ^A 2 ' 

// 18 (0 = 2^(1 + ^A 2 ) 

where f = x - {2/j, A + a 0 ) t. 
Family 11: Ml9 (0 = « 0 -2y^(l + M ) ^ 1+A ^ tanh ^ + ^) 
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K2o(£) = a Q-^V-^{ 1 + ^A 2 ) 



( 



1 + Ay=^coth(^l + 



//19(f) = -2f*(l + ^A 2 ) 



( 



1 + Ay^tanh^ + y^f) 



sec/2 (A + y/-p4) 



( 



1 + Ay^<coth(A + v^f) 



csch(A + y^f) 



) 



where £ = x- {a 0 - 2\[i) t. 
Family 12: w 2i (f) = « 0 + 2y7<csc/z(A + y^/f), 

"22(f) = a 0 Tl2y/fisech(A + yf^C), 

H 21 (Q=fj(l + 2csch 2 (A + ^Z)), 

H 22 (Q = p(l-2sech 2 (A + ^f)), 

where f = x - a 0 t. 

Family 13: u 2 ^(Q = <2 0 +y^(coth(A + y^f) + csch (A + y^Zf ) + Ay = ^), 

"24(f) = a 0 + v / = ^(tanh(A + v^f) + Isech(A + y^f) + Ay^<), 

^23 (f) = +A«(coth(A + y^f) csch(A + y^f) + csch 2 (A + y^f)). 

H 24 (Q = Tp(/tanh(A + y^f) sec/? (A + y^«f) ± sech 2 (A + y/=j%)). 

where f = x - {a 0 - fi A) t. 
Family 14: 1125(C) — a o~V~fi{ coth(A + y/-fiC) ± csch(A + y/-fiC} + Ay^), 

"26(f) = «o~ v/~^( tanh(A + y^f) ±Isech\A + yf^C] + Xy/^ji), 

// 25 (f) = ±f*{coth(A + y/=jtf) csch(A + y^f) ± csch 2 (A + y/=jty), 

// 2 6(f) = ±ft(ltaxih(A + yf^C) sech(A + y^f ) + sech 2 (A + y^f)), 

where f = x - (p A + a 0 ) t. 
Trigonometric Function Solutions for fi > 0: 

Family 15: w 2 7(f) = « 0 + 2 V^( tan (^~V^)~ cot (^~\/^))' 
"28(f) = a 0 +2y/ji(cot(A + y^f)-tan(A + y^f)), 
// 27 (f) = -2^(tan 2 (A->f) + cot 2 (A-^f) + 2), 
// 28 (f) = -2ft (tan 2 (A + y^f) + cot 2 (A + >f) + 2), 

where f = x - dot. 
Family 16: W 29 (f) = « 0 " 2 V^( tan (-A ~ \/^f) ~ cot (A -y^f)), 

"30(f) = « 0 -2 v ^(cot(A + y£f) - tan(A + y^f)), 
// 29 (f) = -2^(tan 2 (A->f) + cot 2 (A -y^f) + 2) , 
//30(f) = -2p (tan 2 (A + >f) + cot 2 (A + y£f) + 2), 
where f = * - a 0 £- 
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Family 17: w 3 i(0 = a 0 +^(ttm(A-^/fit;) ± sec(A-^f)), 

«32(f) = fl 0 +V^(cot(A + v^f) ± csc(A+ v^f)), 

H 3 i(0 = -^(sec^-vX) ± tan(^-^) sec(^l-^)), 

#32(0 = -^(cot(A+ csc(A + ± csc 2 (A+ 

where £ = x - a 0 t. 
Family 18: w 33 (£) = a 0 - ^(tar^A -^f) + sec(A-^77<f)) ' 

m 34 (£) = a 0 - y/fi(cot(A + + csc(A + v^f)), 

#33(0 = -^(sec 2 ^-^) + tan^-y^O sec(A-^)), 

H 3i (Q = -^(csc 2 (A + + cot (A + v^f) csc(A + ^)), 

where <f = x - a 0 t. 
Family 19: m 3 s(0 = « 0 + 2 v / ^( cot (^-\^0 + ^v^)> 

w 3 e(0 = a 0 + 2 V^( tan (A + ) + A-^) , 

H 35 (0 = -2[tcsc 2 (A-^), 

H 36 (Q = -2ftsec 2 (A+^), 

where £ = x- {2/4 X + a 0 ) t. 

(2+2^A 2 ) v ^tan(A- v ^,f) 2 v ^(l+A v ^tan(^l -^)) 

Family 20: w 37 (c) = fln + 7 — \ — H 7 ^ , 

7 i7V ' 0 l+A^Star^A-^) tai^A-^) 

(2 + 2^A 2 ) v /^cot(A + 2^7(1 +A v ^coth(A + v ^Q) 
«38«j-«o+ 1+A ^ coth ( A + ^) + cot(A + ^) 

// 37 (0 = (2A, + a 0 ) L + (2 + W 2 )V-^-y-f) + 2^(l+A^tan(A ^))\ 
V ; V ^ ; ^ l+A^tan^-^) tanh(A-^) J 

1 / (2 + 2^A 2 ) v ^tan(A- v ^0 2^(1 + A ^ tan (A- A 
2\^° 1 + A v /Z7tan(A- v /Z7£) tan^-^) J 

+ 4^ 2 A 2 + 4^-2A^a 0 -2 a o 2 , 

/ (2 + 2^A 2 ) v ^cot(A + v ^) 2yg(l +A v ^cot(A+ 
= {2X, + « 0 ) + ————^—^ + cot( ^ + ^ } J 

1 / (2 + 2^A 2 ) v ^cot(A + ^) 2^7(1 +A v ^cot(A+ ^f)) 



2 ^ 1 + A^7 cot(A + ^Sf) cot (A + ^f) 

+ 4^ 2 A 2 + 4^-2A/i«o-ifl 0 2 , 

where £ = x - {2Xf4 + a 0 K- 
Family 21: u 3 g{Q = a 0 -2 y /fi(cot(A-yffti;) +A V /Z7) : 

w 4 o(0 = a 0 -2 v /^(tan(A + ^0 +A^<), 
H 39 (Q = -2^csc 2 (A- v ^0' 
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H i0 (Q = -2ft sec 2 (A + ^), 
where £ = x - {a 0 - 2X\a) t. 

(2+2/<A 2 ) v ^tan(A- v ^f) 2 v 7Z(l+A v /^tan(,4- v ^)) 

Family 22: u 41 (Q - a 0 - 1+Xy/jitaa ( A _^ tan(^->f) ' 

a (2 + 2ftA 2 ) v ^cot(A + ^.f) 2 v ^(l + A v ^cot(A + y^f)) 

/ (2 + 2^A 2 ) v ^tan(A- v X) 2^(l + A^<tan(A-^)) 
H 41 (0 - (« 0 -2A^) ^ 1+A ^ tan( ^_^ tan(A-^) 

l/ (2 + 2^A 2 ) v ^tan(A- v ^<f) 2y^(l +Xy/fi tan (A-y^)) 
2 I a ° l+Xy^tan(A-y^^) tan^-y^f) 

+ 4^ 2 A 2 + 4(4 + 2X(i a 0 -^a 0 2 , 

H ff) = f 2X)( ( 2 + 2 ^ 2 ) Vficot(A + yjg) 2yg(l + X^2cot(A + ^)) \ 
421 j [a ° l+Ay^cot^ + y^) cot(A + y>f) y 

1 / (2 + 2^A 2 ) v ^<cot(A + v ^f) 2 v ^(l+A v ^cot(A + vK))V 
~ 2 l+I^rotJ^+vK) cotpT^) J 

+ 4^ 2 A 2 + 4f< + 2Xftao--ao 2 , 

where £ = x-(a 0 - 2X(i) i. 
Family 23: u^t;) = « 0 ±2y^sec( J A-y / ^<f). 

«44(f) = a 0 ± 2^/ficsc(A + y/flt;), 

/f«(0 = Kl-2sec 2 (A-v^)), 
J fY 44 (0=^(l-csc 2 (A+ v ^f)), 

where £ = x- dot. 
Family 24: » 4S « ) = a 0 + 2^( 1 + M 2 ) ( ^Lfr-Jjg) J • 

»«<0--^(ttS&)' 

where f = x - {2(i X + a 0 ) t. 
Family 25: MO = «o " ^ + ^ 2 ) ^ 1+Av ^ tan (^) J > 

cot(A + ^) \ 
^l+A v ^cot(A + v /^f) J' 



M 48 (f ) = « 0 - 2^(1 + ^A 
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sec(A-y/ji(;) 



) 



2 



2/4, (1 + ^A 2 ) 



1 + Xy/fitm(A-y/fiO 



esc (A + y/jif;) 



2 



#48 (0 



2p (1 + ^A 2 ) 



l + \y/jicot(A + 



where £ = x - [a 0 - 2Xp) t. 
Family 26: U/&{^) = a 0 ± 2 v /^csc(A- v /^f) . 

«5o(0 = « o ±2 V / ^ SeC 0 4 + \^f)> 

//49(f) =^(l-2csc 2 (^l-^)), 

// 50 («f)=^(l-2sec 2 (A+^)), 

where f = - a 0 £- 
Family 27: m 5 i(£) = « 0 - y/ji (cot (A + csc^-^f) + Ay^<), 

«52(0 = a 0 - V^(ta n (^ + V^O + sec(A + y/jil;) + A V /^), 

= +A<(cot(A-^) csc(A + + csc 2 (A-^<0)' 

H 52 (Q = +1* (tan(A + sec(A + rftf) ± sec 2 (A + ^)), 

where <f = x - {a 0 - fi X) t. 
Family 28: m 53 (£) = a 0 - y/jl (cot(A -y/jit;) ± csc(A-y^£) + A^), 

w 54 (f ) = a 0 - yfji ( tan (A + y^f) ± sec(A + + Xy/p) , 

H 53 (Q = ±fi(cot(A -y/ji$) csc(A - ± csc 2 (A 

// 54 (f) = ±p ( tan (A + y/ft) sec(A + + sec 2 (A + ^)), 

where <f = * - A + <2 0 ) 

Remark 

All the obtained solutions have been checked with Maple by putting them back into 
the original equations and found correct. 

Graphical representations 

Some of our obtained solutions are graphically represented in Figures 1, 2, 3, 4, 5, 6, 7, 
and 8. 

Comparisons with (GIG) -expansion method 

Wang et al. (2008) examined exact solutions of the Variant Boussinesq equations by 
using the (G7G) -expansion method and obtained three solutions. On the contrary by 
using the enhanced (676) -expansion method in this article we have obtained fifty six 
solutions. Furthermore, If we set A = 0, a 0 = ft 0 then our solutions H 7 (<f), wn(£) (Family 
5) coincide with the solution Eq. (1.2) obtained by Wang et al. (2008) for A = 0, A 1 = 0 
and if we set A = 0, a 0 = ji 0 then our solutions u s {&, u 12 (0 coincide with the solution 
Eq. (1.2) obtained by Wang et al. (2008) for A = 0, A 2 = 0. Again if we set A = 0 then our 
solutions //7(f)- #11(0 coincide with the solutions Eq. (1.1) obtained by Wang et al. 
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(2008) for A = 0, A 2 = 0 and if we set A = 0 then our solutions H$(£), H 12 (£) coincide 
with the solutions Eq. (1.1) obtained by Wang et al. (2008) for A = 0, A 2 = 0. Corres- 
pondingly, for similar conditions our solutions of Family 19 and Family 21 are coincide 
with the solution Eq. (1.5) and Eq. (1.6) obtained by Wang et al. (2008). Rest of the 
solutions are freshly derived using enhanced (G'/G) -expansion method. 

Conclusions 

Traveling wave solutions to nonlinear evolution equations arising in mathematical 
Physics are of theoretical importance. In this paper, the enhanced (G'/G) -expansion 
method has been successfully applied for obtaining exact traveling wave solutions of 
Variant Boussinesq equations. It has been shown that the enhanced (G'/G) -expansion 
method is quiet capable and well suited for finding exact solutions. The consistency 
of the method gives this method a wider applicability. With the aid of Maple and by 
putting them back into the original equation, we have assured the accuracy of the 
obtained solutions. Finally, it is worthwhile to mention that the method is straightforward 
and concise and it can be applied to other nonlinear evolution equations in engineering 
and the physical sciences. 
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